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Abstract 

The purpose of this paper is to discuss the evolution of modified 
holographic dark energy with variable G in non-flat Kaluza— Klein 
universe. We consider the non-interacting and interacting scenarios of 
the modified holographic dark energy with dark matter and obtain the 
equation of state parameter through logarithmic approach. It turns 
out that the universe remains in different dark energy eras for both 
cases. Further, we study the validity of the generalized second law of 
thermodynamics in this scenario. We also justify that the statefinder 
parameters satisfy the limit of ACDM model. 
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1 Introduction 

The discovery of the accelerating expansion of the universe is a milestone for 
cosmology which has deep implications for the composition of the universe, 
structure formation and its fate. The expansion of the universe shows that it 
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is not slowing down under normal gravity but accelerating due to an unknown 
component termed as dark energy (DE), having a strong negative pressure 
PQ . There are many pieces of evidence for the existence of this component 
of the universe other than the baryonic and non-baryonic dark matter (DM) 
[2J but it has no clear clue about its identity. 

The most convenient explanation for this expansion is the vacuum energy 
that generates sufficient force to push matter apart described by cosmological 
constant [2j. However, there are two alternative ways used extensively in 
order to explain this behavior. The first approach is the work on different 
DE models such as quintessence [3], K-essence [I], phantom [5], quintom 
[5J, tachyon [7], family of Chaplygin gas [SJ, holographic EH] and new 
agegraphic DE [IT]. Among all these models, holographic DE models are 
widely used. They provide the link of the DE density to the cosmic horizon 
[12] and has been tested through various astronomical observations [13J. 

The idea of holographic DE model (HDE) can be extracted from the 
holographic principle, which states that the number of degrees of freedom of 
a physical system should scale with its bounding area rather than its volume 
[Ti] . Cohen et al. JT5] proposed a relation of ultraviolet (UV) A and infrared 
(IR) cutoffs L due to the limit set by forming a black hole in quantum field 
theory. In their point of view, the total energy of the system having size 
L is bounded by the mass of black hole of the same size. Mathematically, 
it can be written as L 3 p^ < LM^, where represents the vacuum energy 

density and M p = (8tvG)~^ is the reduced Planck mass. Thus one can deduce 
holographic DE [9] 

p A = 3m 2 M p 2 L" 2 , 

here constant 3m 2 is used for convenience and G = G(t). 

The variation of Newton gravitational constant G with cosmic time t has 
been considered for discussing the evolution of DE models. This was also 
used for solving the longstanding problems such as the DM problem, the 
controversies of Hubble parameter value and the cosmic coincidence problem 
(references therein [16J). Additionally, a lot of debate is available in literature 
for the choice of IR cutoff: whether it is a Hubble horizon, or particle horizon, 
or future event horizon for flat FRW universe. It was pointed out by Li [10] 
that the future event horizon is the appropriate choice for IR cutoff which 
favors the current observations. 

The modified and multidimensional theories of gravity (including f(-R), 
f(G), f(i?,G) pi], f(T) [18], Brans-Dicke [19], Horava-Lifshitz [20], Kaluza- 
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Klein (KK) [21]) is the second approach in which a phenomenon to modify 
the gravitational sector or increment of dimension takes place. In case of 
higher dimensional theories, KK theory has attracted many people recently 
to discuss the DE puzzle. It exists in two versions: compact (fifth dimension 
is length like and it should be very small) and non-compact forms (fifth 
dimension is mass like) [22J. In addition, on the basis of the N+l dimensional 
mass of the Schwarzschild black hole [23], Gong and Li [21] derived the HDE 
in extra dimensions (called modified holographic dark energy (MHDE)). 

A marvellous work is available which investigates the non-interacting [T0| 
|25| [26] and interacting [27 \ possibilities of HDE with DM in flat and non-flat 
FRW universes. The thermodynamical interpretation of HDE model with 
different IR cutoff was also investigated by many people [2E1 12H] for non-flat 
FRW universe. The proposal of a statefinder {r, s} was given by Sahni el al. 
[30] for characterizing and differentiating various DE models. Alam et al. [3Tj 
proved that these parameters are a useful tool for describing the properties 
of DE models. Some people [32]- [35] have obtained interesting results about 
HDE with the help of these parameters. Also, MHDE has been considered to 
discuss the evolution of the universe [2U EH]- Recently Sharif et al. [23 EE] 
have explored the evolution of interacting MHDE with Hubble horizon and 
event horizon as an IR cutoff in a flat KK universe and also examined the 
validity of generalized second law of thermodynamics (GSLT). 

This paper is devoted to study the evolution of non-interacting and inter- 
acting MHDE in a non-flat KK universe. The sequence of paper is as follows: 
we discuss the evolution of non-interacting and interacting MHDE in non- 
flat universe in the next section. In section 3, the generalized second law of 
thermodynamics is explored. Section 4 contains the discussion of statefinder. 
We summarize our results in section 5. 

2 Modified Holographic Dark Energy and 
Non-flat Kaluza— Klein Universe 

In this section, we evaluate equation of state (EoS) parameter for the non- 
interacting and interacting MHDE with DM in a compact non-flat KK uni- 
verse [39] whose line element is 




ds 2 = dt 2 - a 2 {t)[ 



dr 2 



+ r 2 {dd 2 + sin# 2 ^ 2 ) + (1 - kr 2 )di)\ (1) 



1 — kr 2 
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where a(t) is the cosmic scale factor responsible for the expansion of the 
universe and k = —1, 0, 1 is the spatial curvature which corresponds to open, 
flat and closed universe, respectively. The Einstein field equation in 4 + 1 
dimensions are 

R^-~gtwR = tfr^, (/*, v = o, 1,2,3,4), (2) 

where R^, g^ u , R, T^ u and k represent the Ricci tensor, the metric tensor, 
the Ricci scalar, the energy-momentum tensor and the coupling constant, 
respectively. Also, we assume that the KK universe is filled with perfect 
fluid whose energy-momentum tensor is given by 

T^v = {p + p)u^u u - g^p, (3) 

where p = p\ is the pressure due to DE, p = p A + p m is the DE plus DM 
(dust like) energy densities and is the five velocity satisfying the relation 
u^u^ = 1, respectively. Using Eqs.flSJ) and (J3j), we get the following Einstein 
field equations for non-flat KK universe: 

a 2 
o k 8tvG 

H + 2H + ? = ~r p - (5) 

Here H is the Hubble parameter and dot is the differentiation with respect 
to time. Equation (j3J) can be written in terms of fractional energy densities 
as 

ft m + A = 1 + ft*, (6) 



where 



_ 8nGp m _ 8nGp A _ k 

m ~ 6tf 2 ' A ~ 6# 2 ' fc ~a 2 # 2 ' [ ) 



The MHDE in iV dimensions is given as 

m 2 (N - l^-iL^M 2 



Pa 



2V t 



N-3 



where A^~i and Vns indicate the area of unit iV-sphere and volume of 
the extra-dimensional space, respectively. In case of the KK cosmology, it 
becomes 

3m 2 n 2 L 2 

* = -tort- (8) 
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Here L is defined for non-flat KK universe as 

L = a(t)r(t), (9) 
where r(t) is obtained through the relation 

1 R e 



lo V 1 " kr 2 (t) a(t) 
Further, its integration gives 

r(t) = sinn (y\f\k\u^ = 



u. 




\k 

R e is the the future event horizon which is defined as 

, , f°° da . , f°° dt 

The time derivative of L yields 



a(t) 

L = HL — cosn ( \/\k\u ) (10) 



with 

( cosu, k = +1, 

cosn( y/\k\u) = < 1, k = 0, 

V\ k \ [ coshw, k = -1. 

Using Eqs.flH]) and ( ITU]) , we obtain the evolution of energy density 

Iran 



Pa = Pa 



cosn(\/\k\S) — A G 



(11) 



where = G = G'i? and prime represents derivative with respect 
to lna(t). By differentiating fractional energy densities &>k and using 
Eq.flTDl. it follows that 

Q' A = 2Q A (1 - H ^^ cosn(y/\k\6) - fi£ = -2fi A (l + J^). (12) 

In the following, we investigate EoS parameter for non-interacting and inter- 
acting cases. 
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2.1 Non-Interacting Case 

The equation of continuity for the KK universe is 

p + 4H(p + p) = 0. 

We consider the universe filled with DE and DM which splits it into two 
equations for DM and DE, respectively, as 



p m + AHp m = 0, 
p A + AH{p A + p A ) = 0. 

In order to eliminate the term including H in Eq. (ll2l) . we use Eqs. 
and (TH: 



(13) 
(14) 



2H 



-4 - 2Q k + Q A (6 - 



2mir 



cosn( y/\k\5)) + A G (1 + fi fe - fi A )] . 



(15) 



Inserting this value in Eq. ffl2|) . it follows that 

2m7r 



r2 * — f2/ 



;i-0 A )(6 



cosn{y/\k\S)) + 2tt k - A G (1 + Q k 



2 - tt A (6 



2mTi 



H 2 y/2lh 



cosn( V \k\5)) + 2Vt k - A G (1 + Vt k 



- «a)]. 



(16) 



Now we want to extract EoS parameter w A in terms of redshift parameter 
z. Integration of the conservation equation for DE gives 



p A ~ a 



-4(1+oj a ) 



We follow the procedure of Li [10] and use Taylor expansion of the DE density 
around the present time ao = 1 as 



In p A = In p A + 



o . din p A 



din a 



In a + 



1 d 2 In pa 
2di\na) 2 



(lna) 2 + 
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where p A is the present value of the DE density. The EoS parameter, up to 
second order expansion, becomes 



_ 1 ( dlnp A 



1 / d 2 In pa 



4 \ din a / 8 \d(lna 



Z = UJq + U)\Z, 



here we have assumed the small redshift approximation, i.e., In a 
z) ~ — z, where 



! 1 ( d In p A 



4 \ a ma 
From (0) and (|T3|) . one can get 



1 d 2 In pa 
Wl = 8d(lna) 2 



Pa 



p mo a 4 f2/ 



1 + fifc — Q\ 

Making use of Eqs.JISI), (USD and (JUD, it follows that 



(17) 
ln(l + 

(18) 
(19) 



3 

-77 + 



2 2H 2 ^2tW A 



1 



1 

4 



fi A |fiO|_(i + fiO; 



2H^ A \ni\ i 

4 

m7r 



m 2 7r 2 



G, 



(20) 



m 2 7r 2 



j (_4 _ o° 



+ 3(l + ^) + iA G (l + ^-0 A ))Jl- 2 ^XI^I 



m 2 n 2 



(21) 



where (0) denotes the present time value of the parameters. 

Finally, we obtain o>a by inserting the above equations in (|17[) . In the 
evolution of u\, the MHDE parameter m plays a crucial role and hence we 
plot cja against m as shown in Figure (1). Here we consider the present 
values of fractional energy densities f2 A = 0.73, = 0.01, H Q = 1 [ID] 
and < A G < 0.07 [16]. Plots in turquoise, blue, pink and green curves 
correspond to values of redshift parameter z = 0.1,0.31,0.5,0.9, respec- 
tively. We observe that for z = 0.1, the EoS parameter attains different 
phases of the universe such as phantom (0.04 < m < 0.45), vacuum DE 
(m = 0.45,1.56), quintessence (0.45 < m < 1.56) and then phantom for 
m > 1.56. However, for z = 0.1,0.31,0.5,0.9, the EoS achieves its maximum 
value —0.6, —1.2, —1.25, —1.3, respectively, and remains in the phantom re- 
gion except z = 0.1. 
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Figure 1: Plot of uj a versus m for non-interacting case. 



2.2 Interacting Case 

Here we evaluate EoS parameter in the interacting phenomenon of MHDE 
with DM by using the same procedure as above. In this case, the equation of 
continuity may be converted into two non-conserving equations for DM and 
MHDE, respectively, i.e., 



p m + AHp 

f) A +4:H(p A +p A ) 



T, 



-T, 



(22) 
(23) 



where T denotes the interaction term which can be taken as [23 

T = 4d 2 Hp A . 



The parameter d is a coupling constant and the selection of its square leads 
to the condition of decay from DE to DM. We can find Uq and uj\ by using 
the above procedure as follows: 



2 



mix 



+ i A 



m 2 TC 2 



n° A \n° k \ -(l + n 



m 2 7r 2 



+ 



mn 



G, 



-4 - 0° 



(24) 



+ 3ft A ) + -A G (l + ft°-ft° A +2cWjWl 



2^l^ fc l 
m% 2 



(25) 



Inserting these values in Eq.( fP7|) . we obtain EoS parameter. In this case, u A 
evolves two constant parameters m and d (interacting parameter). Figure 
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Figure 2: Plot of uja versus m for interacting case. 



(2) shows the plot of u\ against m by setting d 2 = 0.1 and the remaining 
parameters are the same as in the non-interacting case. Notice that the 
maximum values of EoS parameter corresponding to the redshift parameter 
becomes smaller than that of the non-interacting case. 



3 Generalized Second Law of Thermodynam- 
ics 

Now we investigate the validity of GSLT for non-interacting and interact- 
ing MHDE with varying G in the non-flat KK universe. According to this 
law, the sum of entropy of matter inside and at the event horizon remains 
always positive with the passage of time [H]. Thermodynamics of black hole 
plays the role of pillar for thermodynamical interpretation of the universe. 
Bekenstein [42J suggested that, in view of the proportionality relation be- 
tween entropy of black hole horizon and horizon area, the sum of black hole 
entropy and the background entropy must be an increasing quantity with 
time. The first law of thermodynamics gives 

TdS = pdV + dE, (26) 

where T, S, E and p represent temperature, entropy, internal energy and 
pressure of the system, respectively. Splitting this law for DE, DM and 
differentiating with respect to time, we obtain 

s a = pA±A, Sm = p=£+^. (27) 



9 



The volume, temperature and entropy of horizon in KK universe become 

7T 2 L 4 „ 1 „ 7T 2 L 3 



V 



T 



S 



H 



2 1 2nU ~ jj 2G 

Also, we require the following thermodynamical quantities: 

Ti 2 L 4 p A n 2 L 4 p m 



(28) 



Equation (jBJ) can be re-written as 



E„ 



1 



— = -i + — + — . 

Pa ^a ^a 

In view of the above equations, we have 



(29) 



(30) 



S-, 



total 



3m 2 7r 4 L 6 
2G 



1 



+ 



mix 



3 (2ft/ 



)lH 2 



(3- A 



^A ^A 



m 2 7T 2 

4Mf 



2if 4 Q A |Qi 



1 - 



m 2 7r 2 



(31) 



where Stota/ is the sum of three entropies. When we substitute the present 
values of &k, Ac uja (with z < 0.75) and m > 0.1 in the above expres- 
sion, it remains non-negative for both interacting and non-interacting cases 
of MHDE i.e., S tota i > 0. 



4 The Statefinder Diagnostic 

Here we explore the behavior of statefinder parameters in the above men- 
tioned scenario. These parameters have geometrical diagnostic due to their 
total dependence on expansion factor. These are defined for a non-flat KK 
universe as [44J 

r = ^' S= 3(g-^)' (32) 

where VL tot = Q\ + Q m = 1 + ftfc and q is the deceleration parameter defined 
as 

« = -Ji- (33) 
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The statefinder parameters are dimensionless and exhibit expansion of the 
universe through higher derivatives of the scale factor. These are a natural 
companion to the deceleration and Hubble parameters. The pair (r, s) defines 
the well-known ACDM model at the fixed point (r, s) = (f^ o t,0). Moreover, 
r can be expressed in terms of the Hubble parameter as 

r = ^- 3 ~3q-2. (34) 
With the help of Eqs. ([3"3"j) and (13"4"|) . one can write 

r = 2q 2 + q-j ! . (35) 

In the non-interacting case, the time derivative of the deceleration pa- 
rameter becomes 

q = H[q{-2 + A G + 2q) - 40 A (1 + 2u A ) + 20 A w A ]. 

Inserting this in Eq. (l35|) . we obtain 

r = (3-A G )[(l + Q k )(l-^A G ) + 2co A n A } + 4u A Q A 

x (1 + 2uj a - ^A G ) - 2Q A u' A , (36) 

(3 - A G )[(1 + n k )(l - |Ag) + 2u A £l A ] + 4u A tt A 

3[(1 + n k )(i - \A G ) + 2u A n A - ^t] 
Au a VL a (2uj a - |Ag) - 2Vt A u' A - Vt tot 

+ 3[(i + n k )(i - \Kg) + 2u A n A - S ' ( ] 



For the interacting case, we differentiate Eq.p]), using ( fill and ( 1151) . and it 
follows that 

— = 8 + $n k + 8u A n A - 2 a g (i + n k + n A u A ) 

+ 8VL A u A (l + u A ) -2Vt A uj' A + 8d 2 VL A . 
Consequently, the corresponding statefinder takes the form 



11 



r (a) non-interacting r (b) interacting 




Figure 3: Plots of r versus s for (a) non- interacting and (b) interacting cases. 

x (1 + 2u A - ^A G ) - 2Q A u' A + 8d 2 Q A , (38) 

(3 - A G )[(1 + n k )(l - |Ag) + 2co A n A ] + 4u A n A 

3[(1 + n k ){l - ±A G ) + 2u A n A - 
(2u A - |Ag) - 2Q A u' A + 8d 2 Q A - Q tot 

3[(l + fi fc )(l-|AG) + 2a; A n A - ' 

We can easily find a single relation of r in terms of s and draw s — r plane as 
shown in Figure (3) for (a) non-interacting and (b) interacting MHDE. Plots 
in pink, blue, turquoise and green colors are drawn at different physically 
acceptable values of m — 0.21,0.61,0.73,0.91 (as already discussed in [38]). 
respectively, for non-interacting as well as interacting cases. Also, we fix 
z = 0.1, d 2 = 0.1 and recover the ACDM model in both cases. 



5 Concluding Remarks 

We have investigated the behavior of EoS parameter, GSLT and statefinder 
for the MHDE (non-interacting and interacting with DM) with variable G 
correction in non-flat KK universe enclosed by future event horizon. Actually, 
the MHDE exhibits the dynamical nature of the vacuum DE through its 
parameter m. We have evaluated the EoS parameter with respect to m for 
different ranges of z = 0.1, 0.31, 0.5, 0.9. It is found that the interacting and 
non-interacting MHDE behave like a quintom model for a comparatively 
smaller value of the redshift parameter, i.e., z = 0.1 with the assumptions 
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of the present values of the other parameters. For other values of z, it 
evolutes the universe in phantom DE era in view of increasing m. Our 
results about evolution of MHDE with varying G shows compatibility with 
the present observations for flat, non-flat FRW [T61I26] and flat KK universes 
[38] enclosed by future event horizon that it can cross the phantom divide. 

Secondly, we have explored that GSLT is satisfied with z > 0.75 for the 
universe describing phantom evolution. Moreover, we have obtained the evo- 
lution of non-interacting and interacting MHDE in the statefinder plane for 
different best fitted values of m, d (for interacting case) and z = 0.1. The 
trajectories of s-r plane have been achieved with respect to different model 
parameters which is started from right to left. Notice that the parameters s 
and r show decreasing and increasing behavior, respectively, with the phan- 
tom evolution of the KK universe. 
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